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Abstract 

In explicit CP violation scenario of the minimal non-minimal supersymmetric 
standard model (MNMSSM), the possibility of a strongly first-order electroweak 
phase transition (EWPT) is investigated at the one-loop level, where the radiative 
corrections from the loops of the top and stop quarks are taken into account. As- 
suming that the stop quark masses are not degenerate, the radiative corrections due 
to the stop quarks give rise to a CP phase, which triggers the scalar-pseudoscalar 
mixing in the Higgs sector of the MNMSSM. The lighter stop quark need not al- 
ways to have a small mass in order to ensure the strongly first-order EWPT. In the 
MNMSSM with explicit CP violation, it is found that the strength of the first-order 
EWPT depends on several factors, such as the lightest neutral Higgs boson mass 
and the nontrivial CP phase arising from stop quark masses. The effects of these 
factors are discussed. 
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I. Introduction 



For theoretical models to explain successfully the observed baryon asymmetry of the 
universe, it is necessary that they should satisfy the conditions which have been suggested 
by Sakharov several decades ago [1]. The three conditions for dynamically generating the 
baryon asymmetry are: the violation of baryon number conservation, the violation of both 
C and CP, and the deviation from thermal equilibrium. 

The baryogenesis via the electroweak phase transition (EWPT) [2-11], which in prin- 
ciple may satisfy the Sakharov condition, has been studied by many authors. As is well 
known, in order to ensure sufficient deviation from thermal equilibrium, the EWPT should 
be strongly first order, since otherwise the baryon asymmetry generated during the elec- 
troweak phase transition would subsequently disappear. In general, the strength of the 
EWPT is measured by comparing the vacuum expectation value (VEV) of the Higgs field 
at the broken-phase state with the critical temperature. The electroweak phase transition 
is said to be strong if the former is larger than the latter. 

The Standard Model (SM) has already been investigated whether it can realize the 
strongly first-order EWPT. It is found, however, that the SM faces severe difficulties 
to accommodate the desired EWPT, because it cannot make the EWPT strong enough 
unless the mass of the SM Higgs boson is below the present experimental lower bound. 
Thus, in the SM, the EWPT is weakly first or higher order for the experimentally allowed 
mass of the Higgs boson [12-17]. Also, the SM cannot produce the CP violation large 
enough to generate the baryon asymmetry by means of the complex phase in the Cabibbo- 
Kobayashi-Maskawa matrix alone [18]. 

Many scenarios for expanding the SM have been studied in order to formulate the 
idea of baryogenesis via the EWPT [19-33]. Embracing the supersymmetry (SUSY) is 
obviously the most plausible possibility to expand the SM. As the simplest version of 
the supersymmetric SM, the minimal supersymmetric standard model (MSSM) is studied 
widely in the context of baryogenesis via the EWPT. In the MSSM, a strongly first-order 
EWPT is only possible in a confined parameter region where a stop quark is lighter than 
the top quark and the other stop quark is as heavy as the SUSY breaking scale (1-2 TeV). 

The possibility of the EWPT is further examined in other versions of the supersym- 
metric SM, by introducing one or several Higgs singlets to the Higgs sector of the MSSM, 
such as the next-to-minimal supersymmetric SM (NMSSM) with a Za-symmetry [24-27], 
the general NMSSM with broken Z3 [28, 29], the minimal non-minimal supersymmetric 
SM (MNMSSM or nMSSM) with tadpole terms [30, 31], and the minimal supersymmetric 
models with an additional U{1)' [32], to cite some of them. These models are different 
from the MSSM in the sense that they do not always need a light stop quark in order to 
produce a strongly first-order EWPT. We also have studied the MNMSSM with a tadpole 
term elsewhere for the possibility of the EWPT [31], with no CP mixing in its Higgs 
sector. We have found that this model has some parameter regions where a sufficiently 
strong first-order EWPT may occur for electroweak baryogenesis, without requiring a 
light stop quark. 

In this paper, we would like to continue the previous study so as to consider the 
case with explicit CP violation in the Higgs sector of the MNMSSM at the one-loop 
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level. It is known that there is no CP mixing between the scalar and pseudoscalar Higgs 
fields in the MNMSSM at the tree level. At the one-loop level, the MNMSSM may 
have the CP violation in an explicit way, generated by the radiative contributions due 
to the top and stop quark loops. In this explicit CP violation scenario, we examine if 
it is possible to achieve the EWPT that is first order and strong enough to realize the 
necessary electroweak baryogenesis. 

We find that a strongly first-order EWPT is indeed possible in the MNMSSM with 
explicit CP violation at the one-loop level. The trilinear term in the tree-level Higgs 
potential plays an important role, later at the one-loop level, in achieving a strongly 
first order EWPT in the MNMSSM with explicit CP violation. The nontrivial CP phase 
arising from the stop quark masses alters the strength of a first-order EWPT to be even 
stronger a little. The strength of the first-order EWPT depends also on the lightest 
neutral Higgs boson mass and the stop quark masses, but not in a definite way. As the 
lightest Higgs boson mass increases, the strength of the first-order EWPT may either 
increase or decrease. Also, as the two stop quark masses increase, it may either increase 
or decrease. 



II. THE CP VIOLATING HIGGS POTENTIAL 

The Higgs sector of the MNMSSM consists of two Higgs doublet superfields Hi = {H^, H~), 
H2 = {H~^,H2) and a Higgs singlet superfield [34-36]. Keeping only the Yukawa cou- 
plings for the third generation of quarks, the superpotential of MNMSSM may be written 

as 

W = htQH2t% + hQHib% + \NH^eH2 , (1) 

where hq{q = t, b) are the quark Yukawa couplings, A is a dimensionless coupling constant, 
e is the antisymmetric 2x2 matrix with 612 = 1. The left-handed quark superfields 
qL^q = t,b) form the SU(2) doublet = {ti^bi), and the charge conjugate of the 
right-handed qR^q = t, b) quark superfields form the SU(2) singlets and 6^. 

The tree-level Higgs potential, Vq, of the MNMSSM may be decomposed into i^-terms, 
F-terms, the soft terms, and the tadpole terms as 

Vo = Vd + Vf + Vs + Vt^ , 

where 

Vp = \\nm\' + \H2nN\' + \H^eH2\'], 

Vt = -ieN + R.c.) , (2) 

with gi and g2 being the U(l) and SU(2) gauge coupling constants, respectively, and a 
being the Pauli matrices. In Vs, Ax is the trilinear soft SUSY breaking parameter with 
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mass dimension, and m^/^, mj^^, and m^v are the soft SUSY breaking masses. In Vr, ^ is 
the tadpole coefficient. 

Without the tadpole terms, the above tree-level Higgs potential would additionally 
possess a discrete symmetry which would then lead to a massless pseudo-Goldstone boson. 
This is the unwanted axion, because the discrete symmetry is a Peccei-Quinn symmetry. 
Generally, whenever a discrete symmetry is spontaneously broken down, an unwanted 
axion would appear. 

In terms of the physical Higgs fields, the Higgs doublets and the Higgs singlet may be 
expressed as 

Vi + Si + i sin f3Pi \ 
sin/3C+* j ' 

COS/3C+ \ (3) 

{v2 + S2 + ico^(3Pi)e'' ) ' 

X + 53 + zP2 ) e'^ , 

where Si {i = 1-3) are the scalar Higgs fields. Pi {i=l, 2) are the pseudoscalar Higgs fields, 
C"*" is the charged Higgs field, and vi, V2, and x are respectively the VEVs of Si {i = 1-3) 

with V = \Jvi + V2 = 174 GeV and tan/3 = V2/V1. Note that these VEVs are evaluated at 
zero temperature. Thus, they will be denoted in the next section as fi(0), ^2(0) and a;(0) 
to show explicitly the temperature dependence, whereas the VEVs at finite temperature 
will be denoted as ui, V2, and x. 

If complex phases are present in the Higgs sector of the MNMSSM, we would expect 
CP violation in the MNMSSM. The complex phases may explicitly appear from compex 
coefficients and/or from the VEVs of the relevant Higgs fields. In the tree- level Higgs 
potential. Ax-, A, and ^ can be complex. Among them, it can be easily shown that the 
phase of Ax may be absorbed into that of A, without loss of generality. Further, A can 
be made real by adjusting the phases of the Higgs doublets, while ^ can also be made 
real by adjusting the phase of the Higgs singlet. Consequently, the Higgs potential of the 
MNMSSM at the tree level can be made free of any CP violating complex phases. 

The same conclusion can be drawn by applying the tree-level tadpole minimum con- 
ditions upon the tree-level Higgs potential of the MNMSSM. In general, five complex 
phases may arise from the tree-level Higgs potential of the MNMSSM after spontaneous 
breakdown of the electroweak symmetry. Besides the three complex coefficients Axe^"^"^^, 
Ae*"^^, and ^e*'^«, two of the VEVs can also be complex: ^26*^, and xe''^ . 

A pair of tree-level tadpole minimum conditions are derived from the first derivatives 
of the tree-level Higgs potential with respect the two pseudoscalar Higgs fields as 

= 2\vxAx sin 9i , 

= Xv'^Ax sin 2/? sin 9i + 2^^ sin 62 , (4) 

where the two phases 9i and 62 are defined as 9i = (f)x + (pA^ + ^ + 5 and 62 = 30^ + 6. It is 
straightforward that these two equations are satisfied only when both 6*1 = and 62 = 0. 

Therefore, either way, there is no complex phase that can induce CP mixing between 
the scalar and pseudoscalar Higgs bosons in the the tree-level Higgs potential of the 



Hi 



H2 
N 



4 



MNMSSM. Thus, in order to accommodate any complex phases so as to induce the CP 
mixing between the scalar and pseudoscalar Higgs bosons in the present model, one has to 
consider higher-order corrections. The radiative corrections due to the top and stop quarks 
is known to affect significantly the tree-level Higgs sector of supersymmetric models. 

The one-loop effective potential including the radiative corrections due to top and stop 
quarks is given by [37] 




(5) 



where ti {i = 1 ,2) are stop quarks, A is the renormalization scale in the modified minimal 
subtraction scheme, and M.t and M-i. respectively are the top and stop quark masses 
given as functions of the Higgs fields. 

After the spontaneous breakdown of the electroweak symmetry, the top quark mass is 
given by rrit = /iff 2 and the stop quark masses are given by 



rrin 



m. 



=F htJ Afv2 + X'^vfx'^ + 2XAtViV2X cos( 



(6) 



where is the soft SUSY breaking mass and At is the trilinear SUSY breaking parameter 
with mass dimension. 

One may notice that a complex phase (pt is already present in the stop quark masses. 
It is given by (pt = 4>At + 0a + 6' -|- 5, where each phase is respectively originated from 
Ate^'t'^t^ Xe'^\ V2e^\ and xe^^ . 

The CP-odd tadpole minimum conditions with respect to the pseudoscalar Higgs fields 
are modified at the one-loop level as 








3/^,2 

Ax sin Bx ttt^ A sin 0t/(m| , m?J , 



lg^2 ------ t,: 

v4AAf^sin2^sin^i + 2^^ sin 6*9 



Ath? sin 2^ sin (ptfim^-, , "^|) 



(7) 



where the CP violating phase 9i, 62, and (pt are the same as in the tree level, and the 
scale-dependent function /(m^,m^) is defined by 



ml, ml) 
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{ml - ml) 



2i "^l 2i "^l 



A2 



+ 1 



Substituting the first tadpole minimum condition into the second one, one can easily 
see that the second tadpole minimum condition is satisfied by sin ^2 = 0. Further, the 
first tadpole minimum condition alone shows that 61 at the one-loop level is not zero but 
dependent on the other parameters. Consequently, the MNMSSM may eventually have 
one physical CP phase (pt at the one-loop level even if it is free of any complex phase at 
the tree level. We will rename 61 at the one-loop level as (po hereafter, in order to avoid 
any confusion with the tree-level 61 that is zero. 

Now, let us calculate the masses of the five neutral Higgs bosons. They are given 
as the eigenvalues of the symmetric mass matrix M that is calculated from the second 
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derivatives with respect to Si, S2, S3, Pi and P2- In the basis of {Si, S2, Pi, S3, P2), the 
matrix elements of M are obtained as 

+ {mz cos + "^A sin^ /3, 
+ {mz sin + cos^ /3, 

^2 £3 

+ -^m^ sin^ 2p + — , 

y2 f3 

+ -T^^A sin^ 2/5 + — , 

Ax^ X 

I ^o\2„,2 _2 _2\„- /o„„„A5 



— m 4 sin /3 cos P + 2t;A x cos P 

X 
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Moo 




Mil 
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^55 
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00 


M12 
-^^-^ io 
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Mi4 
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14 


-'►^15 
M23 




mL 

Mi, 


M24 






M25 
M34 




Mi, 
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M35 




Mk 


M45 




Ml, 



1) 

—m\ sin P cos^ P + 2t;A^x sin P 

X 



^ .2 



H — sin /3 cos P , 

X 

(9) 

where m| = (5f^+(7|)t>^/2 is the squared mass of the neutral gauge boson, m\ is introduced 
for convenience as 



2 _ AxAacos0o 3m,^AAxcosv.t , ^2 n ..^^ 
^ " sin/3cos/3 ^ 167rVsin3/3cos/3^^ ' ^ ^ 



and Mlj are the radiative corrections due to the top and stop quarks at the one-loop level. 
Explicitly, M*- are given as 



Mil 



Sm^A^x^Al g{m\, m?J 



Stt^w ^ sin /5 (m?_^ — rn^^y ' 

Sm^A^Al ^(m|, m|) Sm^AA^^ log(m|/m? 



t ^ 

STT^f ^ sin^ P (m?_^ — m?^)^ Att'^v'^ sin^ /5 (m?_^ — m?J 

g^2^2 gjj^S ^ ^ j^4 

^ _ Sm^AVA^sin^^t ^(m? , mf^ 



M; 



87rVsin^/3 (m| - m? 



87r2 tan2 /3 (m? - m? ^2 ' 



t2 h' 



,2 ^2 



M; 



SmfX'^A^ sin^ (f)t 9 {^1, J 



Stt^ tan^ P (m| — m| ^ 
^ SmtXxAtAi^Aj^ ^(m|, m|) ^ 3m^AxA,-^ log(m|/mg 
""^^ Svr^t^^ sin^ /5 (m?_^ — m?^)^ Svr^w^sin^/? (m?_^ — m?^ 



^ ^ _ 3mfX^x^AtAi^sm(j)t ^Kv m|) 

Sn^v^sin^P (m| -m|)2 ' 
^ 3m^A^xAg^ ^(m|, m|) 3mgA^xcot/3 ^ , 

_ 3m^A2xAA-iSin0t (?(m|, m|) 



STT^wsin/Jtan/^ (m|— m|)2 
^ Sm^AxA^ A^-^ sin 0^ , m| ) Sm^XxAt sin 0^ log(m| /m| ) 

^ 3m^AAtA,-,A,-, m|) ^ ' 3m^AA,-^ log(m|/m|) 

Svr^f sin/5tan/3 (m?_^ — m?^)2 STr^f sin /9 tan /3 (m?^— m?^} ' 



M 



M 



_ 3m^AA,2A,-^sin0f m|) _ 3m^AAt sin 0^ log(m| /m| ) 

STT^f sin /? tan /5 (m|— m?^)^ Stt^v sin /? tan /? (m|^ — 
,^ _ 3m^A2a;AAt^ sin0t ^(m|, 



87r2tisin^/5tan/5 (m|— m|)2 
3m^A^a;A^ sin^ 0t m|) 
Svr^f sin^/^tan/? (m| — m|)2 



^35 - ^-i— -2 



,2 j^2 



, ^ 3m^A%A,-^sin0, 

87r2tan2/5 (m? - mH^ ' ^ ^ 

where 

Aj^ = At cos 0f + Ax cot P , 

A^-^ = A + Ax cot /3 cos 0t , (12) 
and g{ml,my) is another scale-independent function that is defined as 

c/(m^,m) = ^ + (13) 

m^. — rriy 

The matrix elements of M that are responsible for the CP mixing between scalar and 
pseudoscalar Higgs fields in the (5*1, 5*2, Pi, 5*3, P2) basis are M13, M23, M15, M25, M34, 
and M45. These matrix elements would be zero unless the radiative corrections M* due 
to the top and stop quarks at the one- loop level are taken into account. Moreover, those 
matrix elements of M* for the CP mixing contain one CP phase (pf Thus, the magnitude 
of the scalar-pseudoscalar mixing is proportion to sin0(. There would be no CP mixing 
between the scalar and pseudoscalar Higgs fields in the Higgs sector of the present model 
in the cases of 0^ = 0, and the mixing would be maximal when sin0f = 1. 

The five physical neutral Higgs bosons hi {i = 1-5) would be defined as the eigenstates 
of the mass matrix M. The eigenvalues of M , denoted as m|, [i = 1-5), are their squared 
masses. We sort these five neutral Higgs bosons in the increasing order of their masses such 
that is the smallest eigenvalue and hi is the lightest neutral Higgs boson. They are in 
general given as the mixtures of 5*1, S2, S3, Pi, and P2- If 0t = 0, these five neutral Higgs 
bosons may be classified into three scalar and two pseudoscalar Higgs bosons. Unless (pt 
vanishes, hi {i = 1-5) would not have definite CP parities. 
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In order to measure the size of the CP violation, which arises from the CP mixing 
between the scalar and pseudoscalar Higgs fields, it is preferable to define the following 
parameter, 

p = h{ol,ol,ol,ol,olfi' , 

where Oij =1-5) are the matrix elements of the orthogonal matrix O that diagonalizes 
M. The parameter p can vary between and 1, since Oij satisfy the orthogonality 
condition of J2^=i Ofj = 1. If p = 0, there would be no CP violation in the Higgs sector 
of the NMSSM. On the other hand, if p = 1, CP symmetry would be maximally violated. 
The maximal CP violation that leads to p = 1 takes place when Ofi = = Of^ = 
01 = 01, = 1/5. 



III. THE THERMAL POTENTIAL 

Up to now, we have studied the Higgs sector of the MNMSSM at zero temperature. The 
relevant quantities such as VEVs and various parameters are all defined and derived at 
zero temperature. In order to study the possibility of the EWPT in the present model, 
we now consider the temperature dependence of the model. In particular, the thermal 
effects due to the top and stop quarks are taken into account for the thermal potential 
at finite temperature at the one-loop level. We employ the effective potential method to 
obtain the radiatively corrected thermal potential at the one-loop level. 

At finite temperature T, the radiative corrections due to the top and stop quarks is 
effectively given as [38] 



Vt{^, T)= J2 ^ dxx^ log 1± exp (-Jx^ + mf{^)/T 



(14) 



where ip denotes collectively the Higgs fields, the negative sign is for the stop quarks 
ti{i = 1, 2) and the positive sign for top quark t, and ni (/ = t, t2, ^2) are the degrees of 
freedom for top and stop quarks counting their charge, color, and spin factors. Specifically, 
nt = —12 for top quark and n^. = 6 (z = 1, 2) for stop quarks. The above radiative 
corrections should be added to the potential at the one-loop level at zero temperature 
to yield the finite-temperature one-loop effective potential. The above potential may 
have several minima, of which some may be degenerate. The global minima of the above 
potential is defined as the vacua at T. 

Let us denote the VEVs of the neutral Higgs fields at T as fi(T), thiT), and x{T). 
Hereafter, we will omit the temperature dependence of these VEVs, except the zero- 
temperature values, namely, fi(0), ^2(0), and x(0). In terms of the VEVs at T, the 
finite-temperature one-loop effective potential may be written as 

{V{vi,V2.x,T)) = {Vo) + {Vi) + {Vt) . (15) 
Explicitly, we have 

nr\ 9l ~^ 92 / 2 2\2 I \2/ 22, 22, 2 2\ , 22, 22, 22 

(Vo) = — -{v^ — v^y + X{v^V2 + v^x + v^x) + m^v^ + m2V2 + m^x 
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2\AxViV2X cos 00 ~ 2^^x 




6T 



{Vt) = - j dx log 1 + exp i-^x'^ + m1{v2)/T^\ 

+ — Y J dx x"^ log 1 — exp -y/x^ + m? (t>i, ^2, a;)/?"^] (16) 



where 0o is a redefined phase from 6i at the one-loop level that arises from the CP odd 
tadpole minimum conditions, and the soft SUSY breaking masses at the one-loop level 
are given by 



,2 

2 



mi = - ^cos2l3 - X\x{Oy + v{Oysm^ P) + XAxx{0)tan/3cos(j)o 



Y^(AV + XAtX cot p cos (f)t)fiml^,mf^] 



ml = ^cos2/?- A^(x(0)2 + t;(0)^cos^/3) + AAAx(O)cot/?cos0o 

+ ^(A^ + AA.xtan/5cos0.)/K, -I) " ^ log (^) 

167r2 ) 327r2^ ""^^ ^ 1^ m? J ' 

A ^ 
ml = -A=^i;(O)2 + ^-^t;(OfAAsin2/3cos0o + ^ 

3/z?Af (0)^ cot /3 , , , , ^ , ^ , N 9 9 X / _N 

+ ig^2^(o) (Ax(O)cos/3 + ^sin/3cos0,)/(m|,m|) . (17) 

As mentioned before, fi(0), ^2(0) and x(0) are the vacuum expectation values evaluated 
at zero temperature. 

Notice that there is a trilinear term with Ax in (Vq), which will be shown to enhance 
the strength of the first-order EWPT. Another trilinear term may also be found in (Vr), 
arising from the thermal effects due to the stop quarks, which can be seen by expanding 
the thermal potential using for example the high-temperature approximation. However, 
we would not use the high-temperature approximation but perform the exact numerical 
integration in the thermal potential. 

Now, we inspect the above finite-temperature one-loop effective potential in order 
to study at what temperature it may have two degenerate vacua. This temperature is 
important for the EWPT, because it is the critical temperature, T^, at which the EWPT 
may take place from the symmetric phase to the broken phase. Since vacua are the 
minima of the above potential, the first derivatives of the potential with respect to the 
VEVs should vanish, that is, dV / dip\^=(^^) = 0. 

The first derivative of the above potential with respect to x yields a minimum condi- 
tion, which may be written as 

= 2m3X - 2\AxViV2 cos 00 + 2X^{vl + vl)x - 2^^ 
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— {Xvix + AtV2 COS (pt)f{ml^ , mf^) 
3T^ 2h'fXvi 



X 




x^ + ml/T^)} 



+ 



TTT^^ 2Ti^^^^+ ^tV2 COS (pt) / 

^12 II' 



X 



(18) 



^x2 + m|/T2{l - exp(-y^x2 + m|/T2)} 



In principle, this minimum condition can be used to express x in terms of other vari- 
ables. In practice, however, the minimum condition cannot be solved analytically, since 
it is nonlinear. We use the bisection method to solve it and obtain x. Substituting x into 
{V{vi, V2, X, T)), we can eliminate x and obtain numerically the values of {V{vi, V2, T)). 

By varying the temperature and searching over the (wi, f2)-plane, we investigate 
{V{vi,V2,T)) in order to find out the critical temperature at which it has two degen- 
erate vacua. Let us denote the two degenerate vacua, at the critical temperature, as 
{via,V2a,xa) and {vib,V2b,xb)- The EWPT may take place from the vacuum with sym- 
metric phase to the vacuum with broken phase. Let {via, V2Ai xa) be the symmetric-phase 
state and (f is, t'2B, a^s) be the broken-phase state. The distance between the two vacua, 
at the critical temperature, is given by 



and it is called as the critical VEV [31]. 

The strength of the first-order EWPT, which is driven by a non-equilibrium (out-of- 
equilibrium) condition, depends on the relative size between the critical temperature Tc 
and the critical VEV Vc- In general, the first-order EWPT is said to be strong if Vc > Tc. 
The main aim of our numerical analysis in the next section is therefore to establish the 
possibility of a strongly first-order EWPT in the MNMSSM with explicit CP violation by 
examining Vc and Tc. 



For our numerical analysis, we fix A = 300 GeV for the renormalization scale and = 175 
GeV for top quark mass. The CP symmetry in the MNMSSM at the one-loop level at zero 
temperature is explicitly violated by 0t, arising from the stop quark contributions in the 
effective potential at zero temperature, which is in fact the only source of CP violation. 
Other relevant free parameters besides 0t in the MNMSSM are: tan/3. A, Ax, x{0), ^, m-r. 
At, and T. 

At given temperature, we examine the value of {V{vi,V2,T)) in the (t>i, ^2) -plane, 
where x has been eliminated numerically by means of the minimum condition. This job 




(19) 



IV. NUMERICAL ANALYSIS 
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is repeated by varying the temperature. As an illustration, we plot the equipotential 
contours of {V{vi,V2,T)) in the (fi, f2)-plane in Fig. 1 for T = 156 GeV. The remaining 
parameters are set as 0j = vr/lOO, tan/3 = 10, A = 0.3, Ax = 500 GeV, x(0) = 50 GeV, 
^ = 80 GeV, rriT = 500 GeV, and At = 100 GeV. 

In Fig. 1, one can see that there are three stationary points for {V{vi,V2,T)). The 
values of {V{vi, V2, T)) at these stationary points tell us that the point at around {vi, V2) = 
(170,320) GeV is a saddle point whereas the other two, at around (fi,t'2) = (6,82) GeV 
and (343,472) GeV, are minima. Moreover, the values of {V{vi,V2,T)) at these two 
minima are found to be equal. Therefore, they are degenerate minima, with a potential 
barrier between them. The first-order EWPT may take place from one vacuum to the 
other one through thermal tunnelling. 

The shape of {V{vi,V2,T)) in Fig. 1 is actually a typical example for the first-order 
EWPT. The temperature of T = 156 GeV is the critical temperature Tc that allows the 
phase transition from the symmetric-phase state to the broken-phase state. It is found 
that the minimum point at around {vi, V2) = (6, 82) GeV is the symmetric-phase vacuum 
whereas the minimum point at around (t'i,f2) = (343,472) GeV is the broken-phase 
vacuum. 

In order to estimate the strength of the first-order EWPT for the potential in Fig. 1, 
we need to know Vc- We obtain x = 34 GeV and x = 519 GeV respectively for (wi,f2) = 
(6,82) GeV and (wi,f2) = (343,472) GeV. Thus, the distance between {via,V2a,xa) = 
(6,82,34) GeV and {vib,V2b,xb) = (343,472,519) GeV, where the subscripts A and 
B denote respectively the symmetric-phase state and the broken-phase one, is roughly 
Vc = 708 GeV. Comparing this value with the critical temperature, Tc = 156 GeV, we 
obtain Vc/Tc ~ 4.53. This ratio implies that the first-order EWPT shown in Fig. 1 is 
strong. 

The numerical analysis of Fig. 1 suggests that there is at least one set of free pa- 
rameters of the MNMSSM with explicit CP violation at the one-loop level which allows 
the possibility of a strongly first-order EWPT. For the parameter set of Fig. 1, other 
relevant values are evaluated. The masses of the neutral Higgs boson are calculated to 
be ruh, = 60.91 GeV, mn^ = 94.92 GeV, rrih, = 130.27 GeV, ruh^ = 292.75 GeV, and 
m/jg = 296.03 GeV. The effective size of the CP violation is p = 0.47 x 10~^, which is very 
small, since (pt = vr/100. The masses of the stop quarks are calculated to be m^^ = 512.70 
GeV and m^-^ = 546.24 GeV. They are much heavier than the top quark. 

By examining other sets of parameters, we find that a strongly first-order EWPT is 
possible for a wide parameter space of the MNMSSM with explicit CP violation at the 
one-loop level. After establishing the possibility of the strongly first-order EWPT, we are 
now interested in the dependence of the strength of the first-order EWPT on such factors 
as the CP phase 0t, the mass of the lightest Higgs boson, and the stop quark masses. 

First, in order to study the effect of the CP phase on the strength of the first-order 
EWPT, we take (pt = vr/2, which is maximal since sin0( = 1. The other parameters 
are set for simplicity as Fig. 1: tan/3 = 10, A = 0.3, Ax = 500 GeV, x(0) = 50 GeV, 
^ = 80 GeV, rriT = 500 GeV, and At = 100 GeV. At = 161 GeV, we find that the 
potential has three stationary points with required properties. The result is shown in 
Fig. 2, where the equipotential contours of (V^) are plotted in the (wi, f2)-plane at that 
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temperature. The coordinates in the (wi, f2)-plane of the two degenerate minima of (V^) 
are {via,V2a) = (6,52) GeV with xa = 33 GeV and {vib,V2b) = (356,487) GeV with 
xb = 533 GeV. The strength of the first-order EWPT is Vc/Tc = 4.65. We also have 
p = 0.76 X 10~^. For the masses of relevant particles, we obtain rrih^ = 60.60 GeV, 
ruh^ = 94.92 GeV, = 130.39 GeV, ruh, = 293.38 GeV, and mh, = 296.66 GeV for 
neutral Higgs bosons, and m^^ = 512.95 GeV and m^^ = 546.01 GeV for the stop quarks. 

The numerical results in Fig. 2 are essentially the same as those in Fig. 1, although 
(pt is respectively nearly zero in Fig. 1 and maximal in Fig. 2. The mass spectra of 
neutral Higgs bosons or the stop quarks are almost the same; the critical temperature 
obtained in Fig. 2 is slightly higher than that in Fig. 1; and the strength of the first-order 
EWPT in Fig. 2 is also a little stronger than that in Fig. 1. It is thus hard to draw any 
conclusive relationship between the effect of (pt and the strength of the first-order EWPT, 
by comparing the numerical results of Fig. 2 with those in Fig. 1 alone. We may suggest 
that, in the MNMSSM with explicit CP violation, the CP phase (pt enhances slightly the 
first-order EWPT. This behavior of the present model might be compared with the results 
of other supersymmetric models. In the NMSSM with Z3 symmetry, for example, the CP 
phase (pt reduces a little the strength of the first-order EWPT [27]. 

One might suggest that the strength of the first-order EWPT is dependent on the 
lightest Higgs boson mass. Note that the lightest Higgs boson mass in Fig. 2 is smaller, 
negligibly, than that in Fig. 1. Thus, one can expect that the strength of the EWPT is 
enhanced as the lightest Higgs boson becomes light. This is reasonable in both the SM 
and the MSSM, where the strength of the EWPT may always increase when the lightest 
Higgs boson mass decreases. Also, the non-minimal supersymmetric models with a Higgs 
singlet can accommodate a stronger EWPT for a lighter hi [28,29]. 

In order to clarify whether or not the smaller lightest Higgs boson mass may indeed 
enhance the strength of the first-order EWPT, in the present model, we study by changing 
the parameter values. We take A = 0.2 and ^ = 100 GeV in order to increase the lightest 
Higgs boson mass. The other parameters are the same as Fig. 2. That is, (pt = 7r/2, 
tan/3 = 10, Ax = 500 GeV, a;(0) = 50 GeV, ttit = 500 GeV, and A = 100 GeV. Note in 
particular that the CP phase (pt = vr/2 is the same as Fig. 2. 

At Tc = 191 GeV, we obtain Fig. 3. The potential has two degenerate minima at 
{via,V2a,xa) = (11,11,40) GeV and iviB,V2B,XB) = (686,753,784) GeV. The strength 
of the first-order EWPT is about Vc/T^ = 6.53. Also, we have p = 0.17 x 10~^. The mass 
spectra of relevant particles are: nih^ = 89.96 GeV, = 130.15 GeV, m^g = 144.88 
GeV, m/i^ = 244.18 GeV, and irih^ = 247.14 GeV for the neutral Higgs bosons, and 
m^^ = 512.95 GeV and m^^ = 546.0 GeV for the stop quarks. 

Comparing Fig. 3 with Fig. 2, where nih^ = 89.96 GeV with Vc/Tc = 6.53 in Fig. 3 
while m/i^ = 60.60 GeV with Vc/Tf. = 4.65 in Fig. 2, one might argue that the strength 
of the EWPT is increases as the lightest Higgs boson mass increases. On the other hand, 
by comparing Fig. 2 with Fig. 1, one might obtain an opposite argument. Thus, it is 
likely that the strength of the EWPT in the MNMSSM with explicit CP violation does 
not consistently dependent on the variation of the lightest Higgs boson mass. It can be 
either enhanced or reduced as the lightest Higgs boson mass increases. 

Now, we investigate the dependence of the strength of the EWPT on the stop quark 
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masses. It has been discussed in the hterature that the stop quarks play some crucial role 
in other supersymmetric models within the context of the EWPT. The MSSM always 
need a light stop quark below top quark mass in order to achieve a strongly first-order 
EWPT, whether or not CP violation occurs in its Higgs sector. This behavior of the 
MSSM is understandable, since the thermal contribution of the stop quarks is suppressed 
as they become very heavy. On the other hand, the MNMSSM with explicit CP violation 
is known to have a parameter set, without need of a light stop quark, where the first-order 
EWPT is as much strong as in the case of Fig. 1 in the present model. 

In the present model, there is the trilinear term containing Ax, already at the tree 
level, which leads to a stronger EWPT. As seen in the previous figures, a strongly first- 
order EWPT is allowed without a very light stop quark in the MNMSSM with explicit CP 
violation. We take rriT = 600 GeV in order to increase the stop quark masses. The other 
parameters are the same as Fig. 3: (pt = vr/2, tan/3 = 10, A = 0.2, Ax = 500 GeV, x(0) = 
50 GeV, ^ = 100 GeV, and At = 100 GeV. With these parameter values. Fig. 4 shows 
for Tc = 145 the equipotential contours of {V) in the {vi,V2)-pla.ne with two degenerate 
vacua at {v^a, V2A, xa) = (11, 123, 44) GeV and {vib, f2B, xb) = (641, 708, 748) GeV. The 
strength of the first-order EWPT is about Vc/T^ = 7.66. We obtain p = 0.13 x 10~^, 
and ruh, = 92.66 GeV, rrih^ = 130.15 GeV, rrih., = 145.68 GeV, ruh^ = 244.18 GeV, 
and m/ig = 247.20 GeV for the neutral Higgs boson masses and m^-^ = 610.838 GeV and 
771^2 = 638.84 GeV for the stop quark masses. 

These stop quark masses may be compared with those in Fig. 3, namely, m^^ = 512.95 
GeV and mj^ = 546.0 GeV. The stop quark masses in Fig. 4 are increased but not so 
significantly changed from Fig. 3. Nevertheless, we find that the strength of the first 
order EWPT is noticeably increased. Thus, we may deduce that heavy stop quarks are 
favored by the strongly first-order EWPT, in the MNMSSM with explicit CP violation. 
This behavior distinguishes the present model from the MSSM. 

V. CONCLUSIONS 

We have established the possibility of a strongly first-order EWPT in the MNMSSM with 
explicit CP violation. At the tree level, the Higgs potential of the model can be made 
free from any complex CP phase by suitably rotating the Higgs fields. At the one-loop 
level, as the radiative corrections due to the stop quarks with non-degenerate masses are 
taken into account, the explicit CP violation can be generated by a CP phase (pf We 
have studied the Higgs sector of the model by considering the finite-temperature one-loop 
effective potential, where the thermal effects due to top and stop quarks are numerically 
evaluated by exact integration. We find that some parameter space of the model allows 
a strongly first-order EWPT. 

We also have studied how a strongly first-order EWPT depends on the effective size 
of CP violation, the neutral Higgs boson masses, and the stop quark masses for given 
sets of parameter values. We find that the strength of the EWPT has not a consistent 
dependence on the variation of the mass of the lightest Higgs boson: It can be either 
enhanced or reduced with increasing the lightest Higgs boson mass. We find that the 
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strength of the EWPT increases as the complex CP phase gets large. It also becomes 
stronger as the stop quarks get heavier. 
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Figure Captions 



Fig. 1 : Equipotential contours of (V) at Tc = 156 GeV in the (wi, f2)-plane for 0j = 
7r/100, tan/3 = 10, A = 0.3, Ax = 500 GeV, x(0) = 50 GeV, ^ = 80 GeV, tut = 500 
GeV, and At = 100 GeV. Tlie symmetric-phase vacuum is located at (t'i,f2) = (6,82) 
GeV with x = 34 GeV and the broken-phase vacuum is located at (^1,^2) = (343,472) 
GeV with x = 519 GeV. The strength of the first order EWPT is about Vc/Tc = 4.53. 
We have p = 0.47 x 10-^ and rrih, = 60.91 GeV, ruh^ = 94.92 GeV, ruh^ = 130.27 
GeV, TJifi^ = 292.75 GeV, and nih^ = 296.03 GeV for the neutral Higgs boson masses and 
m^^ = 512.70 GeV and m^-^ = 546.24 GeV for the stop quark masses. 

Fig. 2 : Equipotential contours of {V) at = 135 GeV in the (wi, f2)-plane for (pt = 7r/2, 
tan/? = 10, A = 0.3, Ax = 500 GeV, a;(0) = 50 GeV, ^ = 80 GeV, my = 500 GeV, 
and At = 100 GeV. The symmetric-phase vacuum is located at (t'i,f2) = (6,52) GeV 
with X = 33 GeV and the broken-phase vacuum is located at {vi,V2) = (356,487) GeV 
with X = 533 GeV. The strength of the first order EWPT is about Vc/Tc = 4.65. We 
have p = 0.76 x 10-^ and mh, = 60.60 GeV, = 94.92 GeV, = 130.39 GeV, 
m/j4 = 293.38 GeV, and nihr, = 296.66 GeV for the neutral Higgs boson masses and 
m^^ = 512.95 GeV and m^^ = 546.01 GeV for the stop quark masses. 

Fig. 3 : Equipotential contours of {V) at = 191 GeV in the {vi, f2)-plane for (pt = vr/2, 
tan/? = 10, A = 0.2, Ax = 500 GeV, a;(0) = 50 GeV, ^ = 100 GeV, = 500 GeV, 
and At = 100 GeV. The symmetric-phase vacuum is located at (fi,t'2) = (11, 11) GeV 
with X = 40 GeV and the broken-phase vacuum is located at {vi,V2) = (686,753) GeV 
with X = 784 GeV. The strength of the first order EWPT is about Vc/Tc = 6.53. We 
have p = 0.17 X lO'^, and mh, = 89.96 GeV, = 130.15 GeV, ruh, = 144.88 GeV, 
m/i4 = 244.18 GeV, and rrihr^ = 247.14 GeV for the neutral Higgs boson masses and 
nit^ = 512.95 GeV and m^^ = 546.0 GeV for the stop quark masses. 

Fig. 4 : Equipotential contours of {V) at Tc = 145 GeV in the {vi, W2)-plane for (pt = vr/2, 
tan/3 = 10, A = 0.2, Ax = 500 GeV, x(0) = 50 GeV, ^ = 100 GeV, my = 600 GeV, 
and At = 100 GeV. The symmetric-phase vacuum is located at (fi,f2) = (11, 123) GeV 
with X = AA GeV and the broken-phase vacuum is located at (wi,f2) = (641,708) GeV 
with X = 748 GeV. The strength of the first order EWPT is about Vc/Tc = 7.66. We 
have p = 0.13 X lO'^, and m,,, = 92.66 GeV, m^^ = 130.15 GeV, ruh^ = 145.68 GeV, 
= 244.18 GeV, and rrihr^ = 247.20 GeV for the neutral Higgs boson masses and 
m^^ = 610.83 GeV and m^^ = 638.84 GeV for the stop quark masses. 
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Figure 1: Equipotential contours of (V) at Tc = 156 GeV in the (wi, f2)-plane for 0^ = 
vr/lOO, tan/3 = 10, A = 0.3, Ax = 500 GeV, x(0) = 50 GeV, ^ = 80 GeV, mr = 500 
GeV, and At = 100 GeV. The symmetric-phase vacuum is located at (fi,f2) = (6,82) 
GeV with x = 34 GeV and the broken-phase vacuum is located at (fi,f2) = (343,472) 
GeV with x = 519 GeV. The strength of the first order EWPT is about Vc/Tc = 4.53. 
We have p = 0.47 x 10~\ and rrih, = 60.91 GeV, rrih^ = 94.92 GeV, ruh^ = 130.27 
GeV, rrih^ = 292.75 GeV, and m/^g = 296.03 GeV for the neutral Higgs boson masses and 
m^^ = 512.70 GeV and m^^ = 546.24 GeV for the stop quark masses. 
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Figure 2: Equipotential contours of {V) at Tc = 135 GeV in the (wi, V2)-p\a.ne for 0t = 7r/2, 
tan/3 = 10, A = 0.3, Ax = 500 GeV, x(0) = 50 GeV, ^ = 80 GeV, nir = 500 GeV, 
and At = 100 GeV. The symmetric-phase vacuum is located at {vi,V2) = (6,52) GeV 
with X = 33 GeV and the broken-phase vacuum is located at (t'i,f2) = (356,487) GeV 
with X = 533 GeV. The strength of the first order EWPT is about Vc/Tc = 4.65. We 
have p = 0.76 x 10-^ and tjih, = 60.60 GeV, m^^ = 94.92 GeV, m,,3 = 130.39 GeV, 
m/j4 = 293.38 GeV, and m/^j- = 296.66 GeV for the neutral Higgs boson masses and 
m^^ = 512.95 GeV and m^^ = 546.01 GeV for the stop quark masses. 
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Figure 3: Equipotential contours of {V) at Tc = 191 GeV in the (wi, ?;2)-plane for 0t = 7r/2, 
tan/3 = 10, A = 0.2, Ax = 500 GeV, x(0) = 50 GeV, ^ = 100 GeV, mr = 500 GeV, 
and At = 100 GeV. The symmetric-phase vacuum is located at {vi,V2) = (11,11) GeV 
with X = 40 GeV and the broken-phase vacuum is located at (fi,t'2) = (686,753) GeV 
with X = 784 GeV. The strength of the first order EWPT is about Vc/Tc = 6.53. We 
have p = 0.17 x 10"^, and mn, = 89.96 GeV, = 130.15 GeV, ruh^ = 144.88 GeV, 
m/i4 = 244.18 GeV, and m/^j- = 247.14 GeV for the neutral Higgs boson masses and 
m^^ = 512.95 GeV and m^^ = 546.00 GeV for the stop quark masses. 
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Figure 4: Equipotential contours of {V) at Tc = 145 GeV in the (fi, V2)-p\a.ne for 0t = 7r/2, 
tan/3 = 10, A = 0.2, Ax = 500 GeV, x(0) = 50 GeV, ^ = 100 GeV, mr = 600 GeV, and 
At = 100 GeV. The symmetric-phase vacuum is located at (^1,^2) = (11,123) GeV 
with X = 44 GeV and the broken-phase vacuum is located at (fi,f2) = (641,708) GeV 
with X = 748 GeV. The strength of the first order EWPT is about Vc/Tc = 7.66. We 
have p = 0.13 X 10"^, and mn, = 92.66 GeV, = 130.15 GeV, ruh^ = 145.68 GeV, 
m/j4 = 244.18 GeV, and m/^g = 247.20 GeV for the neutral Higgs boson masses and 
m^^ = 610.83 GeV and m^^ = 638.84 GeV for the stop quark masses. 
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